Overview of the basic process steps of the invention 



Step 10. Load the current CB (cone beam) projection into computer memory. 
Suppose that the mid point of the CB projections currently stored in memory is yis^) . 



t 

Step 20. Finding families of lines for filtering. 

Using either Fig. 15 or, additionally Fig. 17 in case of the improved algorithm, identify 
the required families of lines. Then, choose a discrete set of lines from each family. 



Step 30. Preparing for filtering. 

Parameterize points on the said lines selected in Step 20 by polar angle. Using 
interpolation, compute the derivative of the CB data {d I dq)Dj^{y{q),P) |^^^ at 

points P on the said lines that correspond to discrete values of the polar angle. 



Step 40. Filtering. For each line identified m Step 20 convolve the data for that Ime 
computed in Step 30 with filter Msiny . 



Step 50. Back-projection. For each reconstruction point x back-project the filtered 
data found in Step 40 according to equation (18), Then go to Step 10, unless there are 
no new CB projections to process or image reconstruction at all the required points x 
have been completed. 
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Step 20. Finding families of lines for filtering 



Step 21. From the family of lines J?q choose an equidistant set of lines that 
are parallel to the spiral tangent and that cover the projection of the region of 
interest onto the detector plane located between and F.j (see Fig. 11). 



Step 22. From the family of lines £^ choose a discrete set of lines that are 
tangent to F, and F_| (see Fig. 12). The extreme left point of tangency on Fj 
should coincide with the point where the double tangent line is tangent to 
Fj . Similarly, the extreme right point of tangency on F,, should coincide 
with the point where the double tangent line is tangent to F_, . 



Step 23. From the family of lines £^ choose a discrete set of lines that are 

tangent to and F^j (see Fig. 14, left panel). In both cases the points of 

tangency do not have to extend beyond the projection of the region of interest 
onto the detector plane. In case of the improved algorithm, instead of the lines 
tangent to r2 and r_2 , we choose a discrete (say, equidistant) set of values 

for ^3 on the curves and F_2 and then determine the lines L G by 
solving equations (21), (22). On both curves the points 53 do not have to 
extend beyond the projection of the region of interest onto the detector plane. 
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Step 30 



Step 30. Preparing for filtering 



Step 31 . Fix a line L from the said set of lines 
obtained in Step 20. 






Step 32. Parameterize points on the said line by polar angle / in the plane 
through yis^) and L . 
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Step 33. Choose a discrete set of equidistant values /j that will be used later 
for discrete filtering in Step 40. 
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Step 34. For each /j find the unit vector fij which points from y^s^) 
towards the point on L that corresponds to /j . 
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Step 35. Using the CB projection data Dj^(y(q),&) for a few values of q 
close to Sq find numerically the derivative (5/ dq)Dj{y{q)^&) for all 
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Step 36. Store the computed values of the derivative in computer memory. 
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Step 37. Repeat Steps 3 1-36 for all line 
we will create the processed CB data co 
at y{s^) . 


s L identified in Step 20. This way 
rresponding to the x-ray source located 



Step 40 



Fig. 20 



Step 40. Filtering 



Step 41. Fix a line L from one of the families of lines £^ identified in Step 
20. 
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Step 42. Compute FFT of the values of the said processed CB data computed 
in Step 30 along the said line. 






Step 43. Compute FFT of the filter 1/sin/ 
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Step 44. Multiply FFT of the filter 1 / sin / (the result of Steps 43) and FFT 
of the values of the said processed CB data (the result of Steps 42). 
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Step 45. Take the inverse FFT of the result of Step 44 






Step 46. Store the result of Step 45 in computer memory 
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Step 47. Repeat Steps 41-46 for ail line 
give the filtered CB data (^o ^ Pj ) ^ 
number from which L was selected, m 


s in the said families of lines. This will 
where m stands for the line family 

= 0,1,2. 



Step 50 



Step 50. Back-projection 



Step 5 1. Fix a reconstruction point X , which represents a point inside the 
patient where it is required to reconstruct the image 



Step 52. If belongs to / {x) , then the said filtered CB data affects the 

image at x and one performs Steps 53-58. If Sq is not inside the 

interval P^^ (x) , then the said fihered CB data is not used for image 
reconstruction at X . In this case go back to Step 51 and choose another 
reconstruction point. 



Step 53. Find the projection x of X onto the detector plane DP{Sq) and the 
unit vector /^(^q,^), which points from. >'(^o) towards JC. 



Step 54. Using Figs. 1 1, 13, and the right panel of Fig. 14 identify the lines 
from the said families of lines and points on the said lines that are close to the 
said projection xAfxis above or below > *en in the case of the 
improved algorithm use equations (21), (22) and Fig, 16 to find the filtering 
lines close to the said projection JC . This will give a few values of 



Step 55. With interpolation estimate the value of ^ t^jiSQ^fiiSQ^x)) from the 
said values of ^S^^.Pj) for fij close to /3(Sq,x), 
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Step 56 
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Step 50. Back-projection (continuation) 



Step 56. Compute the contribution from the said fihered CB data to the image 
being reconstructed at the point X by muhiplying ^ ^.^Sq, P{Sq,x)) by 

-c^ (5, P , x)) l[An^ I X - y(sQ ) |] . The appropriate backprojection 

coefficient is selected using Figs. 11, 13, and the right panel of Fig. 14 (see 

also Fig. 15 for the summary). If x is above £2'^ or below , then in the case 
of the improved algorithm use Fig. 17 to find the appropriate backprojection 
coefficient . 
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Step 57. Add the said contributions to the image being reconstructed at the 
point X according to a pre-selected scheme (for example, the Trapezoidal 
scheme) for approximate evaluation of the integral in equation (18). 
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Step 58. Go to Step 51 and choose a different reconstruction point x . If all 
reconstruction points x have been processed, go to Step 59. 




YES 
▼ 

The reconstructed 
image is complete 



